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Abstract The Andreev current and the subgap conductance 
in a superconductor/ insulator/ ferromagnet (SIF) structure 
in the presence of a small spin-splitting field show novel in- 
teresting features fl}. For example, the Andreev current at 
zero temperature can be enhanced by a spin-splitting field 
h, smaller than the superconducting gap A , as has been re- 
cently reported by the authors. Also at finite temperatures 
the Andreev current has a peak for values of the spin-splitting 
field close to the superconducting gap, h ps A. Finally, the 
differential subgap conductance at low temperatures show 
a peak at the bias voltage eV = h. In this paper we investi- 
gate the Andreev current and the subgap conductance in SFF 
structures with arbitrary direction of magnetization of the F 
layers. We show that all aforementioned features occur now 
at the value of the "effective field", which is the field acting 
on the Cooper pairs in the multi-domain ferromagnetic re- 
gion, averaged over the decay length of the superconducting 
condensate into a ferromagnet. We also briefly discuss the 
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1 Introduction 

Ferromagnetism and superconductivity are antagonistic to 
each other's orders, however their interplay can be realized 
when the two interactions are spatially separated. In this case 
the coexistence of the two orderings is due to the proximity 
effect 13, 0, J3J. Experimentally this situation can be real- 
ized in superconductor/ ferromagnet (S/F) hybrid structures. 
The role of the Andreev reflection is central to the prox- 
imity effect since it provides the mechanism for converting 
single electron states from a normal (N) or ferromagnetic 
metal to Cooper pairs in the superconducting condensate 
1 5 1, [6]. During the Andreev reflection process the electron 
incoming to the N/S (F/S) interface is reflected as a hole 
and a charge 2e is transferred across the interface. As a re- 
sult a long-range electron-hole coherence is induced into the 
non-superconducting material. The Andreev reflection man- 
ifests itself in the subgap conductance, i.e. the conductance 
for voltages smaller than the superconducting gap A . In dif- 
fusive N/S systems the subgap conductance shows the zero 
bias anomaly peak due to the impurity confinement and the 
electron-hole interference at the Fermi level [7|, |8|, |9|. 

At a S/F interface the mechanism of Andreev reflection 
is modified compared to the N/S hybrid structures since the 
incoming electron and reflected hole belong to different spin 
bands [ 10 1. Thus, one expects a suppression of the Andreev 
(subgap) current by increasing the exchange field h, which 
is a measure of the spin-splitting at the Fermi level. Recently 
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it was shown by the authors that this intuitive picture does 
not hold always |Jj . If the voltage exceeds some critical 
value the Andreev current of a tunnel ferromagnet/ insula- 
tor/ superconductor (FIS) structure is enhanced by a small 
exchange field h < A reaching a maximum at h ss eV at zero 
temperature. If one keeps the voltage low but now increase 
the temperature, the Andreev current (as well as the full cur- 
rent at this temperature) shows a peak at h s» A. All these 
novel features were exhaustively discussed in [ 1 1. Finally it 
was shown that the subgap conductance of a FIS junction 
at low temperatures and small exchange fields h < A has a 
peak at eV = h Q], El - Thus its measurement can be used 
to determine the strength of a weak exchange or Zeeman- 
like field in the hybrid structure. The latter can be not only 
the intrinsic exchange field of a ferromagnetic alloy |[T2|| but 
also a spin-splitting field created in the normal metal by a 
magnetic induction B (in which case h = HbB, where fis is 
the Bohr magneton) or by a proximity to the ferromagnetic 
insulator material I PI. 

All these predictions were made in [ill for the mono- 
domain FIS hybrid system. The purpose of this work is to 
consider hybrid structure with a multi-domain ferromagnetic 
metal. We present a quantitative analysis of the electron trans- 
port in FIS tunnel structures where a ferromagnetic layer 
consist of two magnetic domains with arbitrary direction 
of magnetization (so called "superconducting spin-valve" 
03)). We show that in this case the aforementioned features 
of the Andreev current and subgap conductance occur at the 
value of the "effective field", which is the field acting on 
the Cooper pairs in the multi-domain ferromagnetic region, 
averaged over the decay length of the superconducting con- 
densate into a ferromagnet ifTSl . 



2 Model and basic equations 

The model of a SF1F2N junction we are going to study is de- 
picted in Fig. [T]and consists of a ferromagnetic bilayer F1F2 
of thickness In = h + h connected to a superconductor (S) 
and a normal (N) reservoirs along the x direction. We con- 
sider the diffusive limit, i.e the elastic scattering length t is 
much smaller than the decay length of the superconducting 
condensate into a ferromagnet — y/W/2h and the super- 
conducting coherence length ^ = ^/W/2A, where 3) is the 
diffusion coefficient and h is the value of the exchange field 
(we set h = kg = 1 and for simplicity we assume the same 
@ in the whole structure). We also assume that the F1F2 and 
F2N interfaces are transparent, while the SFi is a tunnel bar- 
rier. Thus, the two ferromagnetic layers are kept at the same 
potential as the voltage-biased normal reservoir. The F1F2 
bilayer can either model a two domain ferromagnet or an 
artificial hybrid magnetic structure. 

The magnetization of the Fi layer is along the z direc- 
tion, while the magnetization of the F2 layer forms an angle 




Fig. 1 The SFj F2N junction. The interface at x = corresponds to the 
insulating barrier (thick black line). Interfaces at x = l\ and x = In are 
fully transparent, a is the angle between the magnetization directions 
of F] and F2. 

a with the one of the layer Fi. Both magnetization vectors 
lie in the yz plane. Correspondingly the exchange field vec- 
tor in the Fi is given by h = (0,0, h), and in the F2 layer by 
h = (0, h sin a, h cos a), where the angle a takes values from 
(parallel configuration) to it (antiparallel configuration). 

Under these conditions, the microscopic calculation of 
the electric current through the structure requires solution 
of the quasiclassical equation for the 8x8 Keldysh-Green 
function G in the Keldysh x Nambu x spin space in the 
FiF 2 bilayer Q6), ifTTl . 

r G R G K 
G A 



i@d x J=\H,G], G 2 = l, G 



(1) 



Here H = T Z (E — ha) is the Hamiltonian, J = Gd x G is the 
matrix spectral current, a = [a x , <J y , O z ) are the Pauli matri- 
ces in spin space and T z in Nambu space. The R, A and K 
indices stand for the retarded, advanced and Keldysh com- 
ponents (we use the symbols for 8x8 and T for 4x4 ma- 
trices). In Eq. (Q} we neglect the inelastic collision term, as- 
suming 1 12 to be smaller than the inelastic relaxation length 

ma. 

In the Fi region h<7 = ha z and the equation Eq. (Q]) has 
the form 

i$dj = [t z (E — a z h) , G] . (2) 

In the F2 region h<7 = h<J z exp(— io x <x) and it is convenient 
to introduce Green's functions rotated in spin space |fl9l , 

G = U^GU, U =exp(/o; r a/2). (3) 

The rotated function G is then determined by Eq. (0. 

The Eq. (0 should be complemented by boundary con- 
ditions at the interfaces. As mentioned above, we assume 
that the F1F2 and F2N interfaces are transparent and there- 
fore the boundary conditions at x = l\ , 1 12 read 

(4) 
(5) 



\x=h-0 \x=k+0> 



drG\ , „ = d x G\ 

X \X=l[— "* 



G , „ = TV. 



(6) 
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At x = 0, the SFi interface is a tunnel barrier, where the 
boundary conditions are given by the relation |20|, 



Here Gs is the Green function of a bulk BCS superconductor 
defined as 



Gs = i 7 u + x x v, 



(u,v) = (E, iA)/e, e = J(E + iri) 2 -A\ 



(8a) 
(8b) 



W <C 1 is the diffusive transparency parameter [21], W = 
4 /2gflR, and Tj is the Dynes parameter 11221 . In our calcu- 
lations we set small r\ = 10~ 3 4o where Aq is the supercon- 
ducting gap at zero temperature. Below we omit tj in ana- 
lytical expressions for simplicity. 

The electric current through the structure is given by the 
following expression l23ll . Il24ll . 



/ = — / Tr x z J K dE , 
8e Jo 



where J K = {Gd x G) K = G R d x G K + G K d x G A . By neglecting 
non-equilibrium effects, the Keldysh component of Green's 
function is related to the retarded and advanced ones by 



(9) 



G K = G R h -nG A , n = n + + X z ti- , 

1 / , E + eV , E-eV 
n+ = — tanh ± tanh 

2 V 2T 2T 



(10a) 
(10b) 



where n± and T are correspondingly the equilibrium quasi- 
particle distribution functions and the temperature. Below 
we express the advanced Green functions through the re- 
tarded ones using the general relation G A = — X Z G R ^X Z l23l . 

In particular, we are interested in the Andreev current, 
i.e. the current for voltages smaller than the superconducting 
gap due to Andreev processes at the SFi interface. It is given 
by the expression |25|, |26|, 



(11) 



T A = \ [ A n_(E)M s (E)Ref Q dE. 
eR Jo 



where M S (E) = A0(A - \E\)/VA 2 - E 2 is the condensate 
spectral function, &(x) is the Heaviside step function and 
the function /o is the singlet component of / at x = 0. This 
equation is used throughout the article to determine the An- 
dreev transport. We neglect the contribution to the Andreev 
current due to the partial Andreev reflection at the energies 
above the superconducting gap. In the case of strong enough 
tunnel barrier at x = this contribution leads to negligible 
corrections [26|. 

Because of the low transparency of the tunnel SFi bar- 
rier, the proximity effect is weak and the retarded Green 
function can be linearized (we omit the superscript R), 



where / is the 2x2 anomalous Green function in the spin 
space (|/| <C 1) that obeys the linearized equation, 



(7) &% e ? = 2Ef-{ho,f}, 



(13) 



where {•, •} stands for the anticommutator. The general so- 
lution of this equation has the form 



f(x)=f(x)+fy(x)Oy+f z (x)O z , 



(14) 



where / is the singlet component and f Z y are the triplet com- 
ponents with respectively zero and ±1 projections on the 
spin quantization axis 11271 . 

Solving Eq. ( fT3l in the Fi layer we obtain for the com- 
ponents of Eq. ( TBI ). 

2W 

f±(x) = a±cosh(k±x) + - — (ua± — v) smh(k±x) , (15a) 
k± 

2W 

f y (x) = a y cosh(k y x) + — — ua y sinh(k y x) 7 (15b) 
k y 

where f± — / ± f z , a, are the boundary values of /; at x = 
(i stands for +,—,y) and the characteristic wave vectors are 



2(E T h) 



In the F2 layer the general solution has the form, 

fi(x) = bismh[ki(x-h 2 )] , 



(16) 



(17) 



(12) 



where fi are the components of the rotated Green function, 
Eq. (0. Using the boundary conditions at the F1F2 interface, 
Eqs. (01|5]l we obtain a set of six linear equations for the six 
coefficients a,- and that can be solved straightforwardly. 
In particular we are interested in /o = (a+ +a_)/2 which 
enters the equation for the Andreev current, Eq. (fTTT i. Since 
the analytical expression is cumbersome we do not present 
it here. 



3 Results and discussion 

First we briefly review the novel features of the Andreev cur- 
rent for a mono-domain SIFN structure (a — 0), discussed in 
IH. At zero temperature we observe the enhancement of the 
Andreev current at high enough voltages above some criti- 
cal value, see Fig. [2] (a), solid black line. The Andreev cur- 
rent first increases by increasing h, reaches a maximum at 
h R3 eV , and then decays by further increase of the exchange 
field. The enhancement of the Andreev current is due to the 
competition between two-particle tunneling processes and 
decoherence mechanisms. Sharp suppression of the Andreev 
current at h ps eV occurs when the electron-hole coherence 
length y/ W /2eV is cut off by the decay length of supercon- 
ducting condensate into a ferromagnet, = y/W/2h. 

Another feature of the Andreev current, predicted in [Q] 
is the peak at h w A which can be observed only at finite 
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Fig. 2 The /i-dependence of the ratio I A {h)/I A (0) for /] = % and l 2 = %, W = 0.007, a = (solid black line), a = n/2 (dashed blue line) and 
a = k (dash-dotted red line), (a) eV = 0.84, 7 = 0; (b) = 0.84, T = 0.124; (c) eV = 0.34, 7 = 0.254; (d) eV = 0.84, T = 0.254. 



temperatures. The relative height of this peak increases with 
temperature and voltage, see Figs.[2](c) and (d), solid black 
line. In case of large enough values of V and T, we observe 
both the enhancement of the Andreev current by increasing 
h and the peak at h s» A, see Fig. [2] (b), solid black line. 
The peak can be observed only for high enough tempera- 
tures when the upper limit of the integration in Eq. (fTTT > is 
A (at zero temperature the upper limit is eV < A). Then the 
integrand in Eq. (fTTT) has a "dangerous point" at E = h = A. 
This peak can be observed by measuring the full electric 
current through the junction as the single particle current is 
almost independent on h. Note that for the values of temper- 
ature used in our calculations A w Aq. 

Now let us reconsider these features for the two-domain 
situation in case of a — n/2 (dashed blue lines in Fig. |2]l 
and a = n (dash-dotted red line in Fig. [2]). The thickness 
of the F layers is chosen to be l\ = % and h — 9^, l\ short 
enough for the superconducting condensate penetrates both 
ferromagnetic layers and h long enough for the full develop- 
ment of the proximity effect in F1F2 bilayer (at small values 
of I2 the Andreev current is suppressed by the proximity of 
the normal reservoir at x = 1 12) [28]. 

Firs of all, we see that increasing a the features (peaks at 
h s» eV, A ) smear and their amplitude reduces. For a = n we 
do not see any more the enhancement of the Andreev cur- 
rent. Secondly, we see shift of these peaks to the larger val- 
ues of h, which is explicitly seen for a = n/2. The peak at 
h s» eV is shifted to the right (Fig.[2](a), dashed blue line) as 
well as the peak at h ss A (Fig.@(d), dashed blue line). This 
can be explained as follows. The superconducting conden- 
sate penetrates both ferromagnetic layers and feel the "effec- 



tive exchange field" h acting on the Cooper pairs, averaged 
over the length ^ fT51 . The h(a) is gradually reduced as a 
increases from to %. As before the Andreev current peak 
is at h(a) s» A which in the case of a finite a corresponds to 
larger values of the bare h, therefore we observe shift of the 
Andreev current peak to the right. 

Let us now calculate the subgap differential conductance 
Ga = (Ha/cIV at zero temperature. It is known that for a 
diffusive NIS junction the differential conductance at low 
temperatures has a peak ateV—A and a zero bias anomaly 
(ZBA) peak due to the impurity confinement and the electron- 
hole interference at the Fermi level Q, (8], J5). It occurs 
at zero bias since for V = the electron is perfectly retro- 
reflected as a hole during the Andreev reflection process. 
Thus the electron and the reflected hole interfere along the 
same trajectory and the interference effect strongly enhance 
the subgap conductance at zero bias |8|. 

For the FIS structures with h < A the ZBA peak is now 
shifted to the finite voltage eV — h (T), ifTTI . see Fig. [3] solid 
black lines. This can be described as follows. Upon entering 
of the Cooper pair into the ferromagnetic metal the spin up 
electron in the pair lowers its potential energy by h, while the 
spin down electron raises its potential energy by the same 
amount. In order for each electron to conserve its total en- 
ergy, the spin up electron must increase its kinetic energy, 
while the spin down electron must decrease its kinetic en- 
ergy, to make up for these additional potential energies in 
F |29|. Therefore the electron-hole pair in F has now the 
momentum mismatch, i.e. the electron is not perfectly retro- 
reflected. However, if eV — h there is a possibility for exact 
retro-reflection (and interference along the trajectory) of an 
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Fig. 4 The bias voltage dependence of the differential conductance at T = for exchange field h = 0.54 for(a)/j =0.5<^ and (b) l\ —2t,,l2 = 9^, 
W = 0.007, a = (solid black line), a = it/2 (dashed blue line) and a = n (dash-dotted red line). Here G A = ARjGa- 



electron to a hole with a same kinetic energy equal to the 
Fermi energy. 

In case of the two-domain ferromagnetic metal we now 
have the ZBA shift to the "effective exchange field" eV = 
h(a). The "effective field" is smaller than the bare h, h(a) < 
h, and therefore we observe the shift of the differential con- 
ductance peak to the left. We can explicitly see this for a = 
n/2 (Fig. [3] dashed blue lines). For a = n the situation is 
more complicated as the effective exchange field is rather 
small in the antiparallel configuration. For l\ = % we ob- 
serve a broad ZBA peak at V = for a — % for both values 
of h = 0.34 and 0.54 (Fig. [3] dash-dotted red lines). For 
h £, the maximum is shifted from the zero bias, see Fig.|4l 
dash-dotted red lines. 

In the end we notice that SIFFN structures with two- 
domain ferromagnetic interlayer show interesting behavior 
of the heat current through the structure 11281 . It is known 
that in NIS tunnel junctions the flow of electric current is 
accompanied by a heat transfer from the normal metal into 
the superconductor [30 1, [31 1. This happens due to the selec- 
tive tunneling of high-energy quasiparticles out of the nor- 
mal metal in presence of the superconducting energy gap A . 
The heat transfer through NIS junctions can be used for the 
realization of a microcooler, and the important problem is 
to overcome possible limitations of its cooling performance. 



Some of the limitations arise from the fact that nonequilib- 
rium quasiparticles injected into the superconducting elec- 
trode accumulate near the tunneling interface [32|, [33 1. This 
problem can be solved by imposing a local thermal equilib- 
rium in the superconductor electrode by means of a "quasi- 
particle trap", i.e. an additional normal metal layer cover- 
ing the superconductor electrode |34|, Il35l . Another fun- 
damental limitation arise from the Andreev reflection pro- 
cesses: the Andreev current does not transfer heat through 
the N/S interface but rather generates the Joule heating I^V 
which fully dissipates in the normal electrode and dominates 
quasiparticle cooling at low temperatures 1261 . [36|. In order 
to reduce this factor it was proposed to add a ferromagnetic 
interlayer in the NIS structure to suppress the Andreev cur- 
rent and enhance the heat current and cooling performance 

From our studies we can conclude that the ferromag- 
netic interlayer with small enough exchange field will rather 
enhance the Andreev current and suppress the heat current 
(cooling power) through the structure. One need ferromag- 
net with an exchange field (bare or effective in multi-domain 
case) larger than the superconductor gap A to suppress the 
Andreev reflection processes and enhance the cooling per- 
formance. The a-dependence of the heat current in SIFFN 
structures was discussed in [28]. 
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4 Conclusions 

To summarize, we have studied the Andreev current and 
the subgap conductance behavior in SIFFN hybrid struc- 
tures with arbitrary direction of magnetization of the F lay- 
ers. We have revisited all novel features predicted recently 
in the mono-domain SIF system in the presence of a small 
spin-splitting field h [ 1|, namely the Andreev current peaks 
at h « eV at T = and at h ss A for high enough temper- 
ature, and the differential conductance peak at eV = h. We 
have shown that in the two-domain case the aforementioned 
features occur at the value of the "effective exchange field" 
h(a) < h, which is the field acting on the Cooper pairs in the 
multi-domain ferromagnetic region, averaged over the decay 
length of the superconducting condensate into a ferromag- 
net, ^/,. Increasing a from to n one gradually reduce the 
effective field h{a). We also briefly discuss the heat trans- 
port and electron cooling in the considered structures. 

Acknowledgements This work was supported by the Spanish Min- 
istry of Economy and Competition under Project FIS201 1-2885 1-C02- 
02, the Basque Government under UPV/EHU Project IT-366- 07, the 
"Topological Quantum Phenomena" (No.22103002) KAKENHI on In- 
novative Areas, a Grant-in- Aid for Scientific Research (No. 22710096) 
from MEXT of Japan, and JSPS Institutional Program for Young Re- 
searcher Overseas Visits. The work of A.O. was supported by the CSIC 
and the European Social Fund under JAE-Predoc program. 



References 

1. A. Ozaeta, A. S. Vasenko, F. W. J. Hekking and F. S. Bergeret: 
Phys. Rev. B 86, 060509(R) (2012) 

2. A. A. Golubov, M. Yu. Kupriyanov, E. Il'ichev: Rev. Mod. Phys. 
76,411 (2004) 

3. A. I. Buzdin: Rev. Mod. Phys. 77, 935 (2005) 

4. F. S. Bergeret, A. F. Volkov, K. B. Efetov: Rev. Mod. Phys. 77, 
1321 (2005) 

5. A. F. Andreev: Zh. Eksp. Teor. Fiz. 46, 1823 (1964); [Sov. Phys. 
JETP 19, 1228 (1964)]; D. Saint-James: J. Phys. (Paris) 25, 899 
(1964) 

6. B. Pannetier, H. Courtois: J. of Low Temp. Phys. 118, 599 (2000) 

7. A. Kastalsky, A. W. Kleinsasser, L. H. Greene, R. Bhat, F. P. Mil- 
liken, and J. P. Harbison: Phys. Rev. Lett. 67, 3026 (1991) 

8. B. J. van Wees, P. de Vries, P. Magnee, and T. M. Klapwijk: Phys. 
Rev. Lett. 69, 510(1992) 

9. F. W. J. Hekking and Yu. V. Nazarov: Phys. Rev. Lett. 71, 1625 
(1993); Phys. Rev. B 49, 6847 (1994) 

10. M. J. M. de Jong and C. W. J. Beenakker: Phys. Rev. Lett. 74, 
1657 (1995) 

11. M. Leadbeater, C. J. Lambert, K. E. Nagaev, R. Raimondi, and 
A. F. Volkov: Phys. Rev. B 59, 12264 (1999) 

12. T. Kontos, M. Aprili, J. Lesueur, X. Grison, and L. Dumoulin: 
Phys. Rev. Lett. 93, 137001 (2004) 

13. A. Cottet: Phys. Rev. Lett. 107, 177001 (2011) 

14. T. Yu. Karminskaya, A. A. Golubov, and M. Yu. Kupriyanov: 
Phys. Rev. B 84, 064531 (201 1) 

15. F. S. Bergeret, A. F. Volkov, and K. B. Efetov: Phys. Rev. Lett. 86, 
3140 (2001) 

16. K. D. Usadel: Phys. Rev. Lett. 25, 507 (1970) 

17. D. A. Ivanov and Ya. V. Fominov: Phys. Rev. B 73 214524 (2006) 



18. K. Yu. Arutyunov, H.-P. Auraneva, and A. S. Vasenko: Phys. Rev. 
B 83, 104509(2011) 

19. FS. Bergeret, A. F. Volkov and K. B. Efetov: Phys. Rev. B 66, 
184403 (2002) 

20. M. Yu. Kuprianov and V. F. Lukichev: Zh. Eksp. Teor. Fiz. 94, 139 
(1988); [Sov. Phys. JETP 67, 1163 (1988)] 

21. E. V. Bezuglyi, A. S. Vasenko, V. S. Shumeiko, and G. Wendin: 
Phys. Rev. B 72, 014501 (2005); E. V. Bezuglyi, A. S. Vasenko, 
E. N. Bratus, V. S. Shumeiko, and G. Wendin: ibid. 73, 220506(R) 
(2006) 

22. J. P. Pekola, V. F. Maisi, S. Kafanov, N. Chekurov, A. Kemppinen, 
Yu. A. Pashkin, O.-P. Saira, M. Mottonen, and J. S. Tsai: Phys. 
Rev. Lett. 105, 026803 (2010) 

23. A. I. Larkin and Yu. N. Ovchinnikov: in Nonequilibrium Super- 
conductivity, edited by D. N. Langenberg and A. I. Larkin (Else- 
vier, Amsterdam, 1986) 

24. W. Belzig, F. K. Wilhelm, C. Bruder, G. Schon, and A. D. Zaikin: 
Superlatt. Microstruct. 25, 1251 (1999) 

25. A. F. Volkov, A. V. Zaitsev, and T. M. Klapwijk: Physica C 210, 
21 (1993) 

26. A. S. Vasenko, E. V. Bezuglyi, H. Courtois, and F. W. J. Hekking: 
Phys. Rev. B 81, 094513 (2010) 

27. F. S. Bergeret, A. F. Volkov, and K. B. Efetov: Phys. Rev. Lett. 86, 
4096 (2001) 

28. A. Ozaeta, A. S. Vasenko, F. W. J. Hekking and F. S. Bergeret: 
Phys. Rev. B 85, 174518 (2012) 

29. E. A. Dernier, G. B. Arnold and M. R. Beasley: Phys. Rev. B 55, 
15174(1997) 

30. F. Giazotto, T. T. Heikkila, A. Luukanen, A. M. Savin, and 
J. P. Pekola: Rev. Mod. Phys. 78, 217 (2006) 

31. J. T. Muhonen, M. Meschke, and J. P. Pekola: Rep. Prog. Phys. 75, 
046501 (2012) 

32. A. S. Vasenko and F. W. J. Hekking: J. Low Temp. Phys. 154, 221 
(2009) 

33. S. Rajauria, H. Courtois, and B. Pannetier: Phys. Rev. B 80, 
214521 (2009) 

34. J. P. Pekola, D. V. Anghel, T. I. Suppula, J. K. Suoknuuti, A. J. 
Manninen, and M. Manninen: Appl. Phys. Lett. 76, 2782 (2000) 

35. D. Golubev and A. Vasenko: in International Workshop on Super- 
conducting Nano-electronics Devices, ed. by J. Pekola, B. Rug- 
giero, and P. Silvestrini (Kluwer Academic, Dordrecht 2002), p. 
165; larXiv:1204.27191 

36. S. Rajauria, P. Gandit, T. Fournier, F. W. J. Hekking, B. Pannetier, 
and H. Courtois: Phys. Rev. Lett. 100, 207002 (2008) 

37. F. Giazotto, F. Taddei, R. Fazio, and F. Beltram: Appl. Phys. Lett. 
80, 3784 (2002) 



